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ABSTRACT 



Context. Magnetohydrodynamic instabilities can be responsible for the formation of structures with various scales in astrophys- 
ical jets. 

Aims. We consider the stability properties of jets containing both the azimuthal and axial field of subthermal strength. A 
magnetic field with complex topology in jets is suggested by theoretical models and is consistent with recent observations. 
Methods. Stability is discussed by means of a linear analysis of the ideal magnetohydrodynamic equations. 
Results. We argue that in azimuthal and axial magnetic fields the jet is always unstable to non-axisymmetric perturbations. 
Stabilization does not occur even if the strengths of these field components are comparable. If the axial field is weaker than 
the azimuthal one, instability occurs for perturbations with any azimuthal wave number m, and the growth rate reaches a 
saturation value for low values of m. If the axial field is stronger than the toroidal one, the instability shows for perturbations 
with relatively high m. 
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1. Introduction 



The magnetic field plays a key role in the formation and 
propagation of astrophysical jets. Polarization measure- 
ments provide important information about the orienta- 
tion and the strength of the magnetic field, which seems 
to be highly organized on a large scale in many jets (see, 
e.g., Cawthorne et al. 1993, Leppanen, Zensus & Diamond 
1995, Gabuzda 1999, Gabuzda et al. 2004, Kharb et al. 
2008). To explain the observational data, simplified mod- 
els of the three-dimensional magnetic field distribution 
have been proposed (see, for example, Laing 1981, 1993, 
Pushkarev et al. 2005, Laing et al. 2006 ). In particular 
the polarization data suggest that the magnetic field has 
a significant transverse component in the main fration of 
the jet volume. In accordance with the "traditional" in- 
terpretation (see, e.g., Laing 1993), this magnetic struc- 
ture can be associated with a series of shocks along the 
jet, which enhances the local transverse field. However, 
Very Long Baseline Interferometry (VLBI) observations 
of BL Lac 1803+784 (Hirabayashi et al. 1998) revealed 
that the transverse magnetic structure can be formed by 
a large-scale toroidal magnetic field, and indeed recent ob- 
servations of a radio jet in the galaxy CGCG 049-033 also 



reveal a predominantly toroidal magnetic field (Bagchi et 
al. 2007, Laing et al. 2006). 

Magnetic structures with a significant toroidal field- 
component are predicted by several theoretical models 
for jet formation (Blandford & Payne 1982, Romanova & 
Lovelace 1992, Koide, Shibata & Kudoh 1998). In par- 
ticular some models of the jet propagation (see, e.g., 
Begelman, Blandford & Rees 1984) propose that the 
toroidal field decays more slowly than the poloidal one 
as a function of the distance from the central object, 
therefore the toroidal field should eventually become 
dominant. Other generation mechanisms based either 
on the combined influence of turbulence and large-scale 
shear (see, e.g. Urpin 2006) or on dynamo action along 
with magneto-centrifugal and reconnection processes (de 
Gouveia Dal Pino, 2005) also predict the occurence of a 
significant toroidal component. A situation where loops 
of the toroidal magnetic field dominate the field distribu- 
tion would produce important dynamical consequences for 
the jet structure, for instance, providing pressure confine- 
ment through magnetic tension forces (Chan & Henriksen 
1980, Eichler 1993). However, it is well known from plasma 
physics (see, e.g., Kadomtzev 1966) that the toroidal 
magnetic field may cause various magnetohydrodynamic 
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(MHD) instabilities even in simple cylindrical configura- 
tions. 

In general, different types of instabilities can occur 
in astrophysical jets. They can be caused by boundary 
effects, shear, stratification and magnetic fields and are 
also likely to be responsible for the observed morpholog- 
ical complexity of jets. For this reason several analytical 
(see, e.g., Ferrari, Trussoni & Zaninetti 1980, Bodo et al. 
1989, 1996, Hanasz, Sol & Sauty 1999; see also Birkinshaw 
1997 for a review) and numerical (Hardee et al. 1992, Appl 
1996, Lucek & Bell 1996, Min 1997, Kudoh, Matsumoto 
& Shibata 1999) works have been devoted to the study of 
the stability properties of jets. A systematic study of the 
Kelvin-Hclmholtz instability in non-relativistic jets with 
the longitudinal magnetic field has been performed by 
Bodo et al. (1989, 1996) including also the effect of rota- 
tion. In particular several unstable modes including a slow 
mode associated with the magnetic field and an inertial 
mode caused by rotation have been found. Shear-driven 
instabilities can also be very important in jets (Urpin, 
2002). Numerical simulations (Aloy et al. 1999a, 1999b) 
indicate that the radial structure of jets may be more com- 
plex with a transition shear layer surrounding the jet core. 
The thickness of the shear layer depends on the distance 
being around 20% of the jet radius near the nozzle but 
broadening near the head where it is of the order of the 
jet radius. This shear layer may have an important con- 
sequence for various properties of jets including their sta- 
bility (see, e.g., Hanasz & Sol 1996, 1998, Aloy & Mimica 
2008). 

In this paper, we consider the instability associated 
with the topology of the magnetic field in jets. Magnetic 
fields generated by hydrodynamic motions can be topo- 
logically non-trivial and thus prone to pressure-driven or 
current-driven instabilities (see, e.g., Longaretti 2008 for 
review). This problem is well studied in the context of 
experiments on controlled thermonuclear fusion (see, e.g., 
Batcman 1978) and has been discussed in relation to jets 
by a number of authors (see, e.g., Eichler 1993, Appl 1996, 
Appl et al. 2000, Lery et al. 2000). A linear stability analy- 
sis of supermagnetosonic jets has been performed by Appl 
et al. (2000) for different magnetic configurations includ- 
ing a force-free field. It has been shown that the current- 
driven instabilities can occur in the axial field for various 
profiles of the azimuthal field even if these components are 
comparable. The non-linear development of the current- 
driven instability has been analysed by Lery et al. (2000), 
who found that this instability can modify the magnetic 
structure of a jet redistributing the current density in the 
inner part. The current-driven instability of relativistic 
jets and plerions in toroidal magnetic fields has been con- 
sidered by Begelman (1998), who obtained that the domi- 
nant instability in this case are the kink (m = 1) and pinch 
(m = 0) modes. The kink mode generally dominates, 
thus destroying the concentric field structure. Apart from 
the current-driven instabilities, the pressure-driven ones 
caused by the thermal pressure gradient can also occur in 
jets. The modes, corresponding to these instabilities, are 



radially localized and characterized by high wavenumbers. 
Likely, these modes are responsible for generating MHD 
turbulence in the jet core (Kersale, Longaretti, & Pclletier 
2000). The pressure-driven instabilities are well studied in 
laboratory pinch configurations as well. Robinson (1971) 
considers the stability of the pinch with a large ratio of the 
gas and magnetic pressure using the energy principle and 
finds that stability is possible if currents are located out- 
side the main plasma column. The criterion of pressure- 
driven instabilities in cylindrical geometry has been de- 
rived by Longaretti (2003). This criterion is related to 
the Alfvenic modes, which usually require less restrict- 
ing conditions for instability and can even be unstable if 
Suydam's criterion of stability is satisfied. If both desta- 
bilizing factors, current and pressure gradient are present, 
then modes of mixed nature can generally show in a jet. 
An example of these unstable modes has been studied by 
Baty & Keppens (2002), who considered the jet to have 
equal density with the surrounding medium, separated by 
a shear layer. The magnetic field in the basic state was as- 
sumed to be supersonic and helical with a vanishing radial 
component. The linear stability analysis mainly addresses 
the behaviour of kink modes with m = ±1 where m is 
the azimuthal wavenumber. In this basic state, the ther- 
mal pressure gradient in combination with the component 
of the current parallel to the magnetic field produces the 
energy of the perturbations that feed the instability. The 
authors calculated the spectrum of these modes and have 
found that they can be unstable. The unstable modes are 
usually oscillatory and their growth rate exhibits a rather 
peculiar double peak dependence on the axial wavevector. 

Note that numerical simulations often support the idea 
that magnetic jets should be unstable and can generate 
turbulent motions caused by instabilities. For instance, 
Lery & Frank (2000) found that the instabilities seen 
in their jet simulations develop with a wavelength and 
growth time that are well matched by a linear stability 
analysis. Relativistic simulations of extragalactic jets by 
Leismann et al. (2005) and Roca-Sogorb et al. (2008) also 
clearly show that the various irregularities that develop in 
the jet flow can be interpreted as MHD instabilities. 

Our study focuses on the stability of jets with a sub- 
thermal magnetic field where the magnetic energy is lower 
than the thermal energy of particles. This kind of mag- 
netic field can be generated, for instance, by the turbulent 
dynamo action or by hydrodynamic motions (stretching) 
in the process of jet propagation. We show that stabil- 
ity properties of jets with a subthermal field can be es- 
sentially different from those with a highly superthermal 
force- free magnetic field as was considered, for example, by 
Appl et al. (2000). In particular, higher azimuthal modes 
with m » 1 can have the highest growth rate at vari- 
ance with the superthermal case. In our jet model, both 
destabilizing factors - the pressure gradient and electric 
current - are presented simultaneously, and the instability 
occurs under the combined influence of both these factors. 
Therefore, the instability considered has a mixed pressure- 
and current-driven nature. 



A. Bonanno, V. Urpin: MHD instabilities jets 



3 



The paper is organized as follows. Section 2 contains 
the derivation of the relevant equations, and Sect. 3 shows 
the numerical results. Conclusions are presented in Sect. 
4. 



2. Basic equations 

We consider a very simplified model assuming that the 
jet is an infinitely long stationary cylindrical outflow and 
its propagation through the ambient medium is modelled 
by a sequence of quasi-equilibrium states. For the sake of 
simplicity, it is often assumed that the velocity of plasma 
within a jet, V, does not depend on coordinates (see, e.g., 
Appl 1996, Appl et al. 2000, Lery et al. 2000), and we will 
adopt the same assumption in our study. Instabilities of 
the magnetic configurations associated to the electric cur- 
rent are basically absolute instabilities, i.e. they grow, but 
do not propagate. We suppose that this is the case also in 
the rest frame of the jet, and unstable perturbations are 
therefore simply advected with the flow at the jet veloc- 
ity (see, e.g., Appl et al. 2000). Therefore, we treat the 
stability properties in the rest frame of the jet. 

We explore the cylindrical coordinates (s, <p, z) with 
the unit vectors (e a , e v , e z ). The magnetic field is as- 
sumed to be axisymmetric with a non-vanishing z- and 
(^-component. The azimuthal field depends on the cylin- 
drical radius alone, B v = B v (s), but the axial magnetic 
field B z is constant in our model. 

We study the behaviour of MHD modes in the in- 
compressible limit. This approximation is well justified 
for modes with the characteristic time-scale longer than 
the period of sound waves and for subsonic motions (see, 
e.g., Landau & Lifshitz 1981). Perturbations of the density 
caused by these motions are small and can be neglected 
in MHD equations. By making use of the incompressible 
limit in a magnetised gas, one can consider slow modes 
with the growth rate (or frequency) lower than the fre- 
quency of fast magnctosonic waves. Since the typical time 
scale of modes under study is of the order of the inverse 
Alfven frequency, our consideration applies if the Alfven 
velocity is lower than the sound speed or, in terms of the 
plasma /3-parameter, if (3 >• 1 (f3 is the ratio of the gas 
and magnetic pressures). Note that the velocity of jet V 
can be much higher than the Alfven and sound speed. In 
the incompressible limit, the MHD equations read 

dv VP 1 

- + ( ..V). = -— + — (VxB)xB, 



V-v = 0, 

— -Vx(vxB) 

V • B = 0. 



(1) 
(2) 
(3) 
(4) 



In the basic state, the gas is assumed to be in hydrostatic 
equilibrium in the radial direction, then 



VP 



1 

47T 



(V x B) x B. 



(5) 



Stability will be studied by making use of a linear per- 
turbative analysis. Because the basic state is stationary 
and axisymmetric, the dependence of perturbations on t, 
ip, and z can be taken in the form exp (at — ik z z — imp) 
where k z is the wavevector in the axial direction and m 
is the azimuthal wavenumber. Small perturbations will be 
indicated by subscript 1, while unperturbed quantities will 
have no subscript. Then, the linearized Eqs. (l)-(4) read 



av i 



VPi 



1 

Airp 



[(V x Bi) x B + (V x B) x Bi], (6) 



V-V! =0, 



aBi - (B ■ V)«i + (wi • V)B = 0, 
V • Bi = 0. 



(7) 
(8) 
(9) 



Eliminating all variables in favour of the radial velocity 
perturbation vi s , we obtain 



d_ 

ds 

2lo b 



1 / 2 2\( dv ls . V U 



k 2 z (a 2 + uj 2 a )v 1s 



k z WB(l — oi) — 



m(l + A) 
s 2 A 2 



1-- 



aX 



(wAz + 2mws) 



rtiLOAz 



s 2 A 2 



Vis 



1+A 



where 



i / j-fi \ 



KB, 



B u 



91ns 



A 




s 2 k 2 , ' 



This equation was first derived by Bonanno & Urpin 
(2008b) in their analysis of the non-axisymmetric stabil- 
ity of stellar magnetic configurations. For axisymmetric 
perturbations (m = 0), Eq. (10) recovers Eq. (11) of the 
paper by Bonanno & Urpin (2008a). 

Instabilities of the magnetic configurations associated 
to the electric current are basically absolute instabilities, 
i.e. they grow but do not propagate. We suppose that this 
is the case also in the rest frame of the jet, and unstable 
perturbations are therefore simply advected with the flow 
at the jet velocity. Equation (10) represents a non-linear 
eigenvalue problem for a, which can be solved once the 
boundary conditions are given. We assume that vi s should 
be finite at the jet axis. As far as the outer boundary is 
concerned, it is somewhat difficult to formulate a plausible 
boundary condition because actually there is no boundary 
between the jet and the ambient medium. Likely, the jet 
is separated from the ambient plasma by the shear layer 
that has a finite thickness. The effect of this thickness on 
the growth rate of instabilities has been studied by Baty 
(2005) who found that the results are only slightly affected 
by a change of the thickness. For the sake of simplicity, 
therefore, we assume that the shear layer is infinitely thin 
(see also Appl et al. 2000). In a supermagnetosonic jet, no 
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signal can propagate from the jet interior to its surround- 
ings, and we can expect that the instabilities behave as if 
the jet is bounded by a rigid conducting wall. Therefore, 
we can mimic the outer boundary by supposing v\ s = at 
the jet radius s = si. Note that we also tried other outer 
boundary conditions (for example, with vu ^ 0), but this 
does not change the results qualitatively. 

We can represent the azimuthal magnetic held as 



B v = B v0 ip(s), 



(11) 



where B^o is the characteristic held strength and ip ~ 1. 
To calculate the growth rate of the instability, it is conve- 
nient to introduce dimensionless quantities 




0.5 



0.0 



0.0 0.2 



0.4 



0.6 



O.f 



1.0 



x= — , q = k z su r= . 

si ojbo 

Then, Eq. (10) reads 



B, 



ipO 



Fig. 1. Dependence of on x for 
p = 1, n = 2 (dashed line), and p 



p = 1, n = 1 (solid line), 
= 2,n = 1. 



d f dvi s ^ 
dx \ dx x 
2q 2 ip(x) 



+ 



dvu 
dx 



where 



q 2 x 2 



x 



dlnA 



dx 



-q Z 1 



me 
qx 2 



(!-«)■ 



q 2 x 2 
2m/ 



Wis 



m z +q x 



4gVV(g) 
a; 2 (r 2 + / 2 ) : 



f = qs + m 



tp{x) 



q 2 x 2 (T 2 + f 2 ) 
m 2 + q 2 x 2 



0, 



(14) 



We thus solve Eq. (13) for a wide range of the parameters 
and different functional dependence VK^)- 

3. Numerical results 

We assume that the function tp(s), determining the de- 
pendence of the azimuthal magnetic held on the radial 
coordinate, is given by 



= x n e~ p ^ x ~^ 



(15) 



where n and p are the parameters. The function ip goes 
to at the jet axis and is equal to 1 at the surface x = 1 
for any n and p. For the purpose of illustration, we plot 
ip(x) for several values of p and n in Fig. 1. Field prohlc 
(15) with n = 1 corresponds to the magnetic configura- 
tions with non-vanishing current at the jet axis, whereas 
this current goes to for n — 2. The total axial cur- 
rent, integrated over the jet section, is the same for both 
cases. Therefore, the distributions with n = 1 and n = 2 
can be the representatives of the jet models with the cur- 
rent density being relatively smoothed and concentrated 
to the outer boundary, respectively. The function (15) has 
the maximum at x. m = n/p and therefore B v reaches the 
maximum at the jet boundary if n > p or inside the jet 
if p > n. Therefore, the cases p = l,n = 1 and 2 repre- 
sent the jet models with the maximum toroidal held at 
the outer boundary, whereas the held has its maximum 



within the jet for p = 2,n = 1. Because the real dis- 
tribution of the magnetic held is unknown in jets, we will 
use Eq.(15) to mimic different possible dependences of B v 
on x. As far as the axial held is concerned, we assume 
that B z is constant. In real jets, however, B z can depend 
(13£m x as well. Generally, this dependence does not qual- 
itatively change the stability properties (see Bonanno & 
Urpin 2008a), but it can provide an additional destabiliz- 
ing effect if B z changes the sign inside the jets (Robinson 
1971). 

In our jet model, the electric current, j = (c/ 4tt) V x B, 
only has a z-component, whereas B has both z- and 
(^-components. Therefore, the current j can be repre- 
sented as a sum of components parallel and perpendic- 
ular to the magnetic held, j|| and j±, respectively. The 
parallel component does not contribute to the Lorentz 
force. However, even if j± = and the magnetic held is 
force-free, this component alone can produce an unstable 
magnetic conhguration. This type of instability is usually 
called current-driven. On the contrary, j± contributes to 
the Lorentz force, which should then be balanced by a 
pressure force in the basic state. The magnetic configura- 
tions with j± 7^ and jy =0 can be unstable as well, 
and this type of instability is called pressure-driven. In 
our model, both destabilizing factors are presented and, 
most likely, the instability is of the mixed nature. To some 
extent, the relative importance of pressure- and current- 
driven effects can be characterized by the ratio j±/j\\ in 
the basic state. This ratio is high if B z <C B v0 (e *C 1) 
and the pressure-driven effects determine stability proper- 
ties. On the contrary, the current-driven effects are more 



important if B z >• B v0 (e >• 1). In Fig. 2 we plot the 
ratio j±/j\\ as a function of the radius for different val- 
ues of parameters determining the magnetic conhguration. 
Evidently the pressure- and current-driven effects gener- 
ally are comparable if e is of the order of 1 . 

Equation (13) together with the given boundary con- 
ditions is a two-point boundary value problem, which can 
be solved by using the "shooting" method (Press et al. 
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Fig. 2. Ratio of perpendicular and parallel to the mag- 
netic field components of the electic current on x for 
e = 0.1, p = l,n = 1 (solid curve), e = 0.5, p = l,n = 2 
(dashed), e = 1, p = l,n = 2 (dot-dashed), and e = 2, 
p = l,n = 2 (dot-dot-dashed). 

1992). In order to solve Eq. (13), we used a fifth-order 
Runge-Kutta integrator embedded in a globally conver- 
gent Newton-Rawson iterator. We checked that the eigen- 
value was always the fundamental mode because the cor- 
responding eigenfunction had no zero except for that at 
the boundaries. 

In Fig. 3 we plot the dependence of the growth rate 
on the axial wavevector for the magnetic configuration 
with non-vanishing current at the axis and a relatively 
weak axial magnetic field, e — B z /B v q = 0.1. As was 
noted by Bonanno & Urpin (2008), an axial field breaks 
the symmetry between positive and negative values of the 
axial wavevector q in Eq. (13) even if B z is small com- 
pared to B v . However, Eq. (13) still contains some de- 
generacy because it is invariant under (m, q) — > (—to, —q) 
or (to, e) — > (—to, —e) transformation. The instability oc- 
curs only for perturbations with q within a narrow range 
that depends on the azimuthal wavenumber to. For exam- 
ple, perturbations with to = 1 and to = 6 are unstable 
if > q > -70 and -60 > q > -200, respectively. Note 
that the unstable perturbations with positive to should 
have negative q and, on the contrary, if to is negative, 
instability occurs only for perturbations with positive q. 

The growth rate has two clear maxima with the high- 
est maximum corresponding to q ~ —m/e. By the order 
of magnitude, the axial wave-vector of the most rapidly 
growing perturbation can be estimated from the condi- 
tion of magnetic resonance 

UA = vfe ( k ' Ba + 7^) " °" (16) 

Indeed, this equation implies oja oc qe + mip{x) m 0. Since 
4>{x) <~ 1 in our model, the condition loa — corresponds 
to 

q~-m/e. (17) 



Fig. 3. Dependence of the normalized growth rate of in- 
stability on the wavevector q for e = 0.1, p = 1 and n = 1. 
Different types of lines correspond to different azimuthal 
wavenumbers: to = (solid), to = 1 (dashed), to = 2 
(dot-and-dashed) , to = 4 (dot-dot-dashed), and to = 6 
(long-dashed) . 

The most rapidly growing modes turn out to be highly 
anisotropic if the axial field is weak compared to the 
toroidal one: their axial wavelength \ z — 2n/k z ~ 2-kes is 
much shorter than the radial and azimuthal lengthscale. 
The growth rate is fairly high and is of the order of the in- 
verse Alfven time scale. The growth rate slowly increases 
with to and perturbations with a shorter azimuthal scale 
grow faster. The axisymmetric mode (to = 0) turns out 
to be the most slowly growing. 

The model with p = l,n = 2 exhibits a similar be- 
haviour of perturbations (see Fig. 4). As mentioned, elec- 
tric currents are more concentrated near the outer bound- 
ary in this model. The values of q that allow the instability 
are smaller in this case and therefore the corresponding 
vertical wavelengths are longer, but the range of unstable 
q is narrower. For the considered values of to, the growth 
rate is lower approximately by a factor 2. The general 
impression is that the configuration with currents concen- 
trated closer to the outer boundary is more stable than 
that with more uniformly distributed currents. This con- 
clusion qualitatively agrees with the result obtained by 
Robinson (1971) from the hydromagnetic energy princi- 
ple. The author considers a pinch configuration with large 
/3 and finds that the configuration is more stable if large 
axial currents flow outside the main plasma column. 

Note that the equilibrium state of both the config- 
urations considered in Fig. 3 and 4 is characterized by 
the negative pressure gradient that is required for the de- 
velopment of instability (Longaretti 2008). Indeed, using 
Eq. (5) and expression (15) for the toroidal field, we obtain 

Evidently dP j ds < everywhere within the range 1 > x > 
if p = 1 and n = 1,2. The sign of the pressure gradient is 
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Fig. 4. Dependence of the normalized growth rate of in- 
stability on the wavevector q for e = 0.1, p = 1 and n = 2. 
Different types of lines correspond to the same azimuthal 
wavenumbers as in Fig. 2. 

important because it determines the destabilizing effect in 
the so called Suydam's criterion. This criterion represents 
a necessary condition for stability (see, e.g., Longaretti 
2003) and reads in our notations 



sB 2 z fldhy g dP 
47r \h ds J ds 



(19) 



where h = sB z /B v is the magnetic shear. For the equi- 
librium configuration with the toroidal field given by 
Eq. (15), this criterion can be rewritten as 



8x 



(1 - n + pxf - V 2 



p >0. 



(20) 



For the chosen parameters (p = 1, n = 1, 2), the necessary 
condition for stability is not satisfied in some fraction of 
the jet volume (for example, near the outer boundary), 
and the corresponding configurations can generally be un- 
stable. 

Figure 5 shows the dependence of the growth rate on 
q for the model with p = 2, n = 1. In this case, the max- 
imum of the azimuthal field is located within the jet vol- 
ume, at x = 0.5. The maximum value is approximately 1.3 
times higher than the boundary one. It appears that the 
instability grows slightly faster in this type of magnetic 
configurations, but all the main qualitative features of the 
instability are unchanged. Because B v is stronger in this 
model, the instability occurs for larger q as follows from 
Eq. (16). Note that we calculate T only for to < 3 because 
of computational problems at larger m. 

In Fig. 6 we plot the growth rate of instability for the 
configuration with the axial field much weaker than the 
toroidal one, e = 0.01. The distribution of the toroidal 
field is given by Eq. (15) with p = 1 and n = 2. The in- 
stability can occur in this case as well, but only perturba- 
tions with very large q can be unstable. Even for the mode 
m = 1, the growth rate reaches its maximum at q w 200. 



Fig. 5. Dependence of the normalized growth rate of in- 
stability on the wavevector q for e = 0.1, p = 2 and 
n = 1. Different types of lines correspond to different to: 
to = (solid), to = 1 (dashed), m=2 (dot-and-dashed) , 
and to = 3 (dot-dot-dashed). 



m=l 




Fig. 6. Growth rate as a function of q for e = 0.01 and 
distribution of the toroidal field (15) with p = 1 and n = 2. 



For modes with larger to, the maximum growth rate corre- 
sponds to substantially higher values of q (very short axial 
wavelengths). Unfortunately, our numerical code does not 
allow us to resolve the modes with a very short axial wave- 
length, and we were able to extend our calculations only 
to q ss 300. Note that that even perturbations with a q 
this high still grow rapidly and the growth rate is of the 
order of the inverse Alfven time. 

Figure 7 shows the growth rate versus q for a higher 
ratio of the axial and toroidal fields (e = 0.5) and for the 
model with p = 1 and n = 2. In this case, the instability 
occurs at substantially lower values of q and, correspond- 
ingly, at longer axial wavelength. The modes with small 
to grow relatively slowly, but the growth rate increases 
with increasing to. The mode with to = 8 has the growth 
rate comparable to the inverse Alfven time, but the modes 
with larger to are growing even faster. In that sense, the 
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Fig. 7. Growth rate as a function of q for e = 0.5, p = 1 
and n — 2. The curves correspond torn = 1 (solid), to = 2 
(short-dashed), m = 3 (dot-and-dashed), m = A (dot-dot- 
dashed) and to = 8 (long-dashed). 



Fig. 8. Comparison of the growth rate for the cases £ = 1 
(solid) and 2 (dashed). The parameters determining the 
distribution of the toroidal field are p — 1 and n = 2. 



dependence of modes on to in the cases e ~ 1 and e < 1 
seems to be qualitatively similar. 

Figure 8 compares the growth rates for the cases e = 1 
and 2. We considered the model with p = 1 and n = 2, 
but the behaviour is qualitatively the same as for the other 
models. The growth rates are generally lower than in the 
previous cases (e < 1), but the instability is still present. 
A suppression of the growth rate is certainly more pro- 
nounced if e — 2. However, even in this case, Y increases 
relatively rapidly with increasing to and it is likely to reach 
a saturation value ~ 1 at sufficiently large to. The values 
of q that correspond to unstable perturbations are lower 
than in Figs. 3-6. Figure 8 illustrates very well that the in- 
stability cannot be completely suppressed even if the axial 
field is sufficiently strong and e > 1. At large e, the in- 
stability is not very efficient (r 2 > but very small) for 
perturbations with not very large to but those with large 
m can still reach a saturation value of Y 2 ~ 0.1. Indeed 
in the specific case of e = 2 we could not find significant 
growth rate for modes with m < 5, but we found that the 
growth rate is of the order of the inverse Alfven crossing 
time for sufficiently high azimuthal wavenumbers. 

To illustrate how velocity perturbations depend on the 
cylindrical radius, in Fig. 9 the eigenfunctions for one 
of the models shown in Fig. 8 (e — 1) are plotted. The 
eigenfunctions correspond to the azimuthal wavenumbers 
m = 1 (solid), 5 (dashed), and 10 (dot-and-dashed). The 
m = 1 mode is the only mode that does not go to zero at 
the jet axis. Typically, the modes reach their maximum 
closer to the outer boundary, and an increase of m results 
in the eigenfunctions that are more and more localized. 
For large to, motions in the core region of the jet are sub- 
stantially suppressed. This sharpness of the eigenfunctions 
is the main reason of the problems in numerical calcula- 
tions for high to. 




0.0 0.2 0.4 0.6 0.8 1.0 



x 

Fig. 9. The eigenfunctions for the model with s = 1, p = 
1, and n = 2. The curves correspomd to to = 1 (solid), 
to = 5 (dashed), and m = 10 (dot-and-dashed). 

4. Discussion 

We considered the instability that can occur in jets if the 
magnetic field has a relatively complex structure with non- 
vanishing azimuthal and axial components. A magnetic 
field with complex topology is suggested by both the the- 
oretical models for jet formation and the mechanisms of 
generation of the magnetic field. Our study considers only 
jets with subthermal magnetic fields because the the sta- 
bility properties of jets with superthermal fields can be 
essentially different (Appl et al. 2000). 

The considered model is very simplified and does not 
take into account many physical factors that can be impor- 
tant in the evolution of real jets. For instance, we model 
the transition from the jet to the external medium in terms 
of the infinitely thin boundary layer. In fact, the bound- 
ary can be somewhat extended and the velocity can change 
smoothly within the transition layer. Shear has a destabil- 
ising influence on jets (Urpin 2002), and the shear-driven 
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instability can often be faster than instabilities associated 
to electric currents. An interplay of the shear-driven ef- 
fects with the magnetic instabilities will be considered 
elsewhere. Another important assumption concerns the 
topology of the background magnetic field. Generally, the 
magnetic geometry of jets can be much more complicated 
than that considered in our paper and can include mag- 
netic structures of various scales. Likely, the influence of 
these structures on the stability of jets can be studied only 
by using 3D numerical modelling. One more essential as- 
sumption of our model concerns the lack of hydrodynam- 
ical motions in the basic state of the jet in a co-moving 
frame. The presence of motions with the velocity com- 
parable to the sound speed (for example, rotation) can 
drastically change the stability properties. 

It turns out that the magnetic jet is always unsta- 
ble if the azimuthal field increases to the outer boundary. 
This does not depend on the ratio of the axial and az- 
imuthal fields. However, the character of the instability 
can be essentially different for different e. If the axial field 
is weaker than the azimuthal one and e < 1, the growth 
rate tends to a saturation value of the order of the inverse 
of a few Alfven crossing times for m > 1. Our code does 
not allow us to calculate eigenfunctions and eigenvalues 
for very large to in this case, but Figs. 3 and 4 indicate 
that the growth rate shows a saturation trend for m^> 1. 
Then the fastest growing mode has the axial wavevector 
k z <~ — m/es, which is determined by the condition of 
magnetic resonance. This wavevector corresponds to the 
axial wavelength A z = 27r/|fc z | <~ (2ne)(s/m). At small e, 
the axial wavelength of the most rapidly growing pertur- 
bations should be very small. Therefore, if the longitudinal 
magnetic field is weaker than the azimuthal one, we ex- 
pect that the instability generates structures in jets with 
a very short length scale in the ^-direction. Note that this 
behaviour can cause problems in numerical modelling of 
the instability because a very high resolution in the axial 
direction is required. 

When the ratio of the axial and azimuthal fields be- 
comes comparable to or greater than 1, the behaviour of 
the instability as a function of to is qualitatively similar, 
although the growth rate is significantly smaller. In this 
case perturbations with relatively small to are substan- 
tially suppressed, but we still find a numerical evidence 
for a saturation trend as we increase to, although it is 
reached for much higher values of to than in the e < 1 
case. It is important to stress that even m ~ 100 in the 
azimuthal direction is still a macroscopic scale, and dissi- 
pative effects are negligible in this case. 

The conlusion that the instability can arise even if e > 
1 with a characteristic azimuthal wavenumber to ^ 1 is at 
variance with the widely accepted opinion that magnetic 
configurations should be stabilised at B z <~ B v . Note that 
the similar conclusion has been obtained by Goedbloed & 
Hagebruk 1972) for the magnetic configuration with the 
constant pitch, B v /sB z = const. 

The dependence of the growth rate on q has generally 
a very particular shape with two maxima. The origin of 



these maxima in T 2 is caused by the resonant nature of 
unstable modes. The resonance is associated with the last 
term on the l.h.s. of Eq.(13) (or Eq.(10)) which is pro- 
portional to / 2 /(r 2 + / 2 ) 2 . For the sake of simplicity, we 
can consider the case p — 0, n = 1. Then, ip/x = 1 and 
the function / (Eq.(14)) characterizes a departure of q 
from the resonant value given by Eq.(17). This resonant 
dependence reaches its maximum value (equal to l/2r 2 ) 
at / 2 = T 2 . At given T, this condition can be fulfilled for 
two values of q, 



These two values of q correspond to two maxima in T 2 (q). 
Note that the gap between the maxima ~ 2r/e increases 
with decreasing e and does not depend on to, which agrees 
qualitatively and quantitatively well with the results of 
our numerical calculations. 

Note that the detailed physical nature and many prop- 
erties of the instability investigated in our paper are not 
well understood yet. Our study also does not allow us to 
answer the question what type of the magnetic configura- 
tion is formed because of the development of instability. To 
answer this question, one needs 3D numerical simulations. 
However, it is possible that the much better understand- 
ing of the nature of the magnetic field in jets is required 
to answer this question. For example, the resulting large 
scale magnetic field can be formed by balancing the rate of 
field decay caused by instability and the rate of generation 
owing to some mechanism (e.g., dynamo). It would be im- 
portant to investigate the consequences of our findings for 
jets by performing more realistic numerical simulations, 
and we hope to address this question in the near future. 
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